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Abstract —In this paper, we characterize the performance of 
relay-assisted underwater wireless optical code division multiple 
access (OCDMA) networks over turbulent channels. In addition 
to scattering and absorption effects of underwater channels, we 
also consider optical turbulence as a log-normal fading coefficient 
in our analysis. To simultaneously and asynchronously share 
medium among many users, we assign a unique optical orthog¬ 
onal code (OOC) to each user in order to actualize OCDMA- 
based underwater network. The most significant challenge in 
underwater optical communication is in the ability to extend 
the short range of its coverage. In order to expand the viable 
communication range, we consider multi-hop transmission to the 
destination. Moreover, we evaluate the performance of a relay- 
assisted point-to-point UWOC system as a special case of the 
proposed relay-assisted OCDMA network. Our numerical results 
indicate significant performance improvement by employing in¬ 
termediate relays, e.g., one can achieve 32 dB improvement in the 
bit error rate (BER) of 10“® using only a dual-hop transmission 
in a 90 m point-to-point clear ocean link. 

Index Terms —OCDMA, serial relaying, chip detect-and- 
forward, log-normal fading, BER performance, turbulent under¬ 
water channel, underwater wireless optical communications. 

1. Introduction 

NDERWATER wireless optical communication (UWOC) 
systems are receiving growing attention for various un¬ 
derwater applications. As opposed to their traditional coun¬ 
terparts, i.e., acoustic communications, they have three main 
superiorities: higher bandwidth, lower latency and better se¬ 
curity. Therefore, UWOC system can be regarded as an 
alternative to meet the requirements of high speed and large 
data underwater communications such as imaging, real-time 
video transmission, high throughput sensor networks, etc 
However, presently, UWOC is suitable only for ranges that are 
typically less than 100 m which hinders its extensive usage. 
This drawback is due to the fact that UWOC suffers from three 
main impairing effects: absorption, scattering and turbulence 
which cause loss, inter-symbol interference (ISI) and fading 
on the received optical signal, respectively. 

Many studies have been focused on characterizing ab¬ 
sorption and scattering effects of different water types, both 
theoretically and experimentally @,0 Based on the experi¬ 
mental results reported in Tang-Dong-Zhang simulated 

UWOC channel by means of Monte Carlo (MC) approach 
considering absorption and scattering effects Q. They also 
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fitted a double gamma function (DGE) to this impulse response 
and evaluated BER numerically, however in the absence of 
turbulence effects. Recently, Akhoundi-Salehi-Tashakori 0 
proposed a cellular OCDMA-based UWOC network where a 
specified optical orthogonal code (OOC) 0 is assigned to 
each user. Eurthermore, potential applications and challenges 
of such a network is elaborated in ©• 

Although many valuable studies have been performed to 
characterize and mitigate turbulence-induced fading in free- 
space optical (ESO) communications 0-(n) its impairing 
effects on the performance of UWOC systems have received 
relatively less attention. However, some relevant and useful 
researches have recently been reported in the literature to 
characterize underwater fading statistics. Eor example in (H 
an accurate power spectrum has been derived for fiuctuations 
of turbulent seawater refractive index. Eurthermore, Rytov 
method has been used in (T^ to evaluate the scintillation index 
of optical plane and spherical waves propagating in underwater 
turbulent medium. And in (E), the on-axis scintillation index 
of a focused Gaussian beam has been formulated in weak 
oceanic turbulence and by considering log-normal distribution 
for intensity fiuctuations the average BER in such systems is 
evaluated. 

It has been shown that the three aforementioned impairing 
factors in underwater medium are incremental functions of 
distance 0, O- This behaviour hinders on the system 
performance for longer ranges. To overcome this limitation 
and therefore to extend the viable communication range, in 
this paper we propose multi-hop transmission over turbulent 
underwater channel. By this scheme, we divide a relatively 
long communication distance with severe absorption, scatter¬ 
ing and fading effects to shorter ones, each with a much 
reduced absorption, scattering and fading effects; and therefore 
with acceptable performance. Eig. 1 depicts an explanatory 
practical example of our proposed topology, namely under¬ 
water optical sensor network. Underwater sensor networks 
will acquire essential roles in investigating climate changes, 
disaster prevention, in monitoring biological, biogeochemical, 
evolutionary and ecological changes in the sea, ocean and 
lake environments, in pollution monitoring, and in helping to 
control and maintain oil production facilities. In this scheme, 
each sensor should ultimately upload its own data to the optical 
base transceiver station (obtsQ Since these sensors are 
battery powered, their low-power signals will not withstand to 

Uhe OBTS is located in the cell’s center, and is connected to the other cells’ 
OBTSs via fiber optic. Detailed experimental prototype of this underwater 
OBTS and the cellular network topology are presented in 0 
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Fig. 1. An illustrative practical example of relay-assisted OCDMA-based 
underwater wireless network: underwater optical sensor network. 


reach a relatively far-located OBTS due to severe absorption, 
scattering and fading effects of the channel. In this circum¬ 
stance, a comparatively simple relay can be implemented 
near this network to detect the received data and forward the 
detected data to the next relay (or to the OBTS in the last hop). 
Further, each sensor transmits its data based on OOC-OCDMA 
technique Q, p5|-fT^ to introduce multiple access capability 
while reducing its interfering effect on the other sensors’ 
signals. In general, in such difficult circumstances the relay- 
assisted strategy is a reasonable solution as opposed to design¬ 
ing an exclusive OBTS nearby the network, since relay nodes 
employment is much simpler and straightforward than OBTSs 
0. Moreover, different cells’ OBTSs are sufficiently distant 
to provide longer range communications and also to reduce 
implementation costs. Hence, to develop temporary repeaters, 
implementation of simple relay nodes with narrow field of 
view (FOV) receivers and narrow/wide beam lasers/LEDs is 
more interesting than complicated omnidirectional OBTSs. 
Also the required processes in relays are much simpler than 
those of OBTSs, which involve different layers of network 
processes 0 0 In addition, our study of relay-assisted 
topology covers the more practical scenario of cooperative 
communication where different users operate as repeaters for 
each other. 

The rest of this paper is organized as follows. In Section 
II, the channel and system models in the context of relay- 
assisted OCDMA network are investigated. In Section III, 
some simplifying assumptions are made in order to approxi¬ 
mate the BER of a relay-assisted underwater OCDMA net¬ 
work. Using Gaussian approximation to calculate the chip 
error rates (CERs) of different hops, the end-to-end BERs for 
both uplink and downlink transmissions are obtained. Further, 
we clarify how to obtain BER performance of a single-user 
relay-assisted point-to-point UWOC link as a special case 
of the aforementioned up-or downlink BERs. In Section IV, 
numerical results for various configurations are presented, and 
finally we conclude the paper in Section V. 

II. Channel and System Model 
A. Channel Model 

Through propagation of optical beam in underwater 
medium, interactions between each photon and seawater par¬ 
ticles consist of absorption and scattering. Energy loss of 
non-scattered light, due to absorption and scattering can be 
characterized by absorption coefficient a (A) and scattering 


coefficient 6(A), respectively. Further, extinction coefficient 
c(A) =a(A) + 6(A) describes the total effect of absorption 
and scattering on energy loss. The values of these coefficients 
can vary with source wavelength A and water types |[T|. It 
has been shown in that absorption and scattering have 
the lowest effect at the interval 400 nm < A < 530 nm, so 
UWOC systems employ the blue/green region of the visible 
light spectrum to transmit and receive data. In Q and |T9| the 
channel impulse response has been simulated based on MC 
approach with respect to the both absorption and scattering 
effects. In this paper the fading-free impulse response between 
any two ith and jth nodes is denoted by To attain this 

impulse response, we simulated the channel similar to jl] and 
|T9| based on MC approach in order to thoroughly include 
absorption and scattering effects. 

In UWOC systems, optical turbulence is the main cause 
of fading on the propagating optical signal through turbulent 
seawater. Optical turbulence occurs due to random variations 
of refractive index. These random variations in underwater 
medium mainly results from fluctuations in temperature and 
salinity | [2Q| . To characterize turbulence effects, we multiply 
ho ij {t) by a (positive) multiplicative fading coefficient hij 
0-0 Weak oceanic turbulence can be modeled with log¬ 
normal distribution HD as; 




2ihij 


27rcr^. 


zcxp 


8cr^ 


( 1 ) 


where (iXij and (j\.. are respectively the mean and variance of 
the Gaussian distributed log-amplitude factor Xij = \hi{hij). 
To ensure that the fading coefficient conserves energy, we 
normalize fading amplitude such that ^^ij] = 1, which 
implies /iXij = ~^Xi Therefore, in describing the fading 
statistics, it demands finmng the dependency of log-amplitude 
variance to the ocean turbulence parameters where we 
will discuss in the remaining of this subsection. 

The scintillation index of a light wave with intensity lij = 
Io,ijhij is defined by |[^, | [I^ ; 


in which lo^ij is the fading-free intensity. Scintillation index 
can be considered as a threshold for separating the weak 
and strong turbulence regimes, i.e., when this parameter is 
smaller than unity the channel is in weak turbulence and vice 
versa (8), (T3), (2T). It has been shown in GD that under the 
assumption of weak turbulence the scintillation index of plane 
and spherical waves can be calculated as; 
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where 0 = 1 and 0 for plane and spherical waves, respectively. 
ko = 27r/A, do and n denote the wave number, propagation 
distance and scalar spatial frequency, respectively. Moreover, 
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Fig. 2. (a) Block diagram of the relay-assisted OCDMA-based underwater 

wireless network; (b) transmitted signal of the nth user in its 0th time slot. 


is the power spectrum of turbulent fluctuations which 
has the form GD, ( 13 ; 


(k) =0.388 X [l + 2.35 (k7?)^/^ 


X ^ - 2u;e“"^^s<5\ ( 4 ) 

where e is the rate of dissipation of turbulent kinetic energy per 
unit mass of fluid, 77 = 10“^ m is the Kolmogorov micro-scale, 
Xt is the rate of dissipation of mean-square temperature and w 
is the relative strength of temperature and salinity fluctuations 
which takes on values in the interval [—5, 0]. Other parameters 
in 0 are as At = 1.863 x 10“^, ^4^ = 1.9 x 10~^. A^tq = 
9.41 X 10"^ and 5 = 8.284(k 7?)'^/^ + 12.978(Kr?)^ (^, 

It can be shown that for turbulent channels with log-normal 
fading distribution the scintillation index is related to the log- 
amplitude variance as cr|.. = exp(4cr^.,) — 1 Therefore, 
for relatively weak oceanic turbulence the scinmlation index 
(which is comprehensively evaluated in |T3| ) directly yields 
the log-amplitude variance and consequently specifles fading 
coefficients’ probability density function (PDF) in ([T]). 

It is worth noting that given a set of conditions, e.g., the 
link range, the water type, inherent optical properties (lOPs), 
etc., each of the aforementioned three impairing effects of the 
channel can be a serious concern. For example, as the water 
turbidity increases the absorption and scattering coefficients 
increase too. Therefore, in turbid harbor waters absorption 
and scattering are very important and should be treated, e.g., 
using multi-hop transmission. On the other hand, the severity 
of turbulence, which is measured by the scintillation index, 
depends on the values of the parameters xt, ^ and w. Based 
on the experimental results in p 2 | the typical values of 5 are 
between 10 “^m^/s^tol 0 “^ ^l^o xt takes the values 

between 7 x 10“^^ K^/s to 10“^ K^/s. Hence, depending on 
the values of these parameters, strong turbulence can occur at 
ranges as short as 10 m and as long as 100 m p3| , which 
impressively differs from atmospheric channels where strong 
turbulence distances are on the order of kilometers. Therefore, 
alleviating such a strong turbulence requires more studies. 


B. System Model 

We consider a collection of M underwater users which 
are located nearby each other and are communicating with 


a relatively distant receiver or OBTS. To simultaneously and 
asynchronously share the medium among the various users, 
a specific OOC code with weight W and length F = Ti^lTc 
is assigned to each user, where and Tc are bit and chip 
duration times, respectively. To guarantee that auto-and cross¬ 
correlation constraints and Ac are equal to unity, maximum 
allowable number of users is upper bounded by; 


M < 


F-l 
W{W -1) 


( 5 ) 


where \_x\ is the integer portion of the real value x 

As stated in the previous section absorption, scattering 
and fading effects are distance-dependent and their severity 
increases for longer ranges. This issue motivates the idea of 
placing intermediate nodes, namely relay nodes, between the 
users and the destination to alleviate the above mentioned 
disturbing effects. Therefore, as illustrated in Fig. 2(a) we have 
located N intermediate relays in a serial topology between the 
users and the destination to increase the viable communication 
range. This topology is also suitable when there is no line-of- 
sight to the destination. In order to simplify relays’ structure, 
chip detect-and-forward (CDF) strategy is employed at the 
relay nodes. In this scheme, each relay simply detects the 
received chips based on an appropriate threshold and forwards 
the detected chips to the next (or final) node. It is obvious 
that chip detection process in the first relay Ri is affected by 
multiple access interference (MAI); however, this interference 
does not affect the chip detection process in the remaining 
relays and the destination, since each of these nodes only 
receives the transmitted signal by the previous node. 

As it is shown in Fig. 2(b), the nth user’s 0th time slot 
transmitted OCDMA-based pulse shape, (t), has the form; 

F-l 

(^)=E (6) 

7=0 

(n) 

where Pc is the nth user’s chip power and is related 
to the average transmitted power per bit of this user as 
Pivl = And = {CF,dy\-,C^F-i) is the 

nth user’s OOC code with length P and weight W p~5| . Pt^ (t) 
is a rectangular pulse with duration Tc and unit amplitude, 
i.e., Ptc{^) = ), where n(t) is a rectangular pulse 

with unit amplitude in the interval [—1/2,1/2]. Therefore, the 
transmitted OCDMA-based signal of the nth user, Sn{t), can 
be expressed as; 

oo 

Sn{t)= E 

k=—oo 

oo F—1 

k= — oo 7=0 

in which G { 0 , 1 } is the nth user’s transmitted bit on its 
kth time slot, representing on-off keying (OOK) modulation. 
Furthermore, the received optical signal from the nth user to 
Pi, after having propagated through the channel with impulse 
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response of represented as; 

VnHt) = Sn{t)*h^nlit) 

oo F—1 

= E T„), (8) 

k= — oo j=0 

where (^) is the aggregated channel impulse response of 
the first hop from the nth user to the first relay, = 

PTc{t) * ^o^ni(^)’ * denotes the convolution operation. 

Tn is the nth user’s delay time with respect to the first relay 

receiver’s reference clock. In this paper we consider to be 
chip synchronous with an integer multiple of Tc eg Note 
that each of the M users produces similar contribution in the 
first relay. Hence, the total received optical signal at the front 
end of the first relay’s detector can be summarized as; 

M 

= (9) 

n=l 

Moreover, various types of noises, i.e., background light, dark 
current and thermal noise all affect the system performance. 
In this paper we adopt Poisson distributed background light 
and dark current and Gaussian distributed thermal noise | [23l . 

We should emphasize that our system model and analytical 
derivations are general and valid for when the transmitters 
use collimated laser beams or diffusive light emitting diodes 
(LEDs). However, from the practical point of view, since the 
users can be mobile, tracking and pointing issues suggest to 
use diffusive links for the first hop. In other words, in order 
to reduce the pointing and tracking difficulties, it is preferable 
to use LEDs or lasers with enough wide divergence angles 
as transmitters of the first hop, during both up-and downlink 
transmissions. Also due to the same reason, receivers with 
larger EOV are preferable for the first hop. On the other hand, 
since relay nodes are fixed, the pointing and tracking troubles 
decrease for the intermediate hops. Therefore, collimated links 
using narrow beam divergence lasers are more suitable for 
these hops to better save energy and reduce the channel tempo¬ 
ral spreading. Eurther, if the destination node is a fixed simple 
receiver, collimated link is more favorable to better save energy 
and mitigate ISI. Einally, for the depicted topology in Eig. 1, 
where destination is an OBTS, due to the physical properties 
of OBTS 0 the transmitter of the last hop must have relatively 
a wide beam divergence angle to better illuminate the separate 
photodetectors of the OBTS. 

HI. BER Performance Analysis 

To evaluate BER, we assume the first user as the desired 
user and the remaining users as undesired interfering users. 
Hence, the average BER of the desired user under the as¬ 
sumption that all users transmit bits “0” and “1” with equal 
probability, can be evaluated as; 

PE = ^Pbe{m + lPbe{0\l), (10) 

in which P 5 e(l| 0 ) and P 5 e( 0 |l) are bit error probabilities 
when bits “0” and “1” are transmitted, respectively. To obtain 
P^e(l|0) and P6e(0|l), we first make the following two 
assumptions: 


1) No Inter-Symbol Interference (ISI): We assume that all 
of the users transmit with a chip rate in which the main portion 
of the energy of the received optical signal r(t) = Ptc(^) * 
ho{t) is confined in the interval (0,Tc), i.e., 

pTc pO poo 

/ r(t)dt:^ / / r{t)dt. (ii) 

do d — oo d Tq 

Under this assumption we can conclude that ISI has a negli¬ 
gible effect on the system performance. Eurther, in order to 
more accurately take into account the loss due to absorption 
and scattering, we simulate the channel impulse response by 
MC approach and we guarantee O by adopting users’ data 
rate not to exceed the predetermined value indicated by O- 
Taking into account that the channel delay spread is much 
smaller than Tc, the channel loss coefficient on a given chip 
can be determined as; 


jy PrMdt 



( 12 ) 


It is worth mentioning that as the link range, water turbidity 
and transmitter beam divergence angle increase the channel 
imposes more attenuation (loss) and time spreading on the 
propagating optical signal. On the other hand, as the receiver 
EOV increases the channel spreading time increases too. In 
these conditions, the assumption in ( pTj ) leads to larger values 
for Tc and hence limits the transmission data rate to smaller 
values. 

2) Identical Chip Power: We assume that all users have 
the same transmit power as Pc. Eurthermore, each relay after 
detecting the received chips, forwards each “ON” chip with 
power Pc. 

In what follows, we apply the aforementioned two as¬ 
sumptions to analyze the BERs of the uplink and downlink 
transmissions. 


A. Uplink BER Analysis 

In this section we focus on the 0th time slot transmitted bit 
of the first (desired) user. All the other users transmit their 
signals asynchronously which cause MAI on the desired user. 
Therefore, with respect to the aforementioned two assumptions 
and ([^, we can express the uplink received signal at Pi as; 

= PjyN^'’ E - jT, - Ti) 

i=o 

M c» F-1 

+ PeE h^nl E (13) 

n=2 k= — oo j=0 

which is composed of two main parts: desired and the un¬ 
desired MAI signals. Let us define a = (o^i, Q( 2 , as 

the interference pattern occurred on the pulsed mark chips of 
the first user’s OOC, where aq is the number of interferences 
on the qth pulsed mark chip of the desired user’s OOC. Eur¬ 
thermore, we define I as the total number of interferences that 
occurred on all the chips of the desired user, i.e., I = 

The random variable (RV) /, which specifies the number of 
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interfering users on the desired user, has a Binomial distribu¬ 
tion function. The joint probability distribution function of the 
interference pattern can be represented as p4| , flSl ; 


Pi{l, a) 


n I'M -iww^Vr 

^ A 2 fA 2Fj 

(14) 


Therefore to evaluate or equivalently Pheibo^lbo^), we 
need to evaluate the BER conditioned on I and a. Then, the 
unconditioned error probability can be evaluated as; 


Pbe{bi^^\bi^^) = 


M-1 

E E 

1=0 de^i 


(15) 


where is the binary complement of the 0th time slot 
transmitted bit of the first user and 0/ is the set of all inter- 
fere nce p atterns with I interfering users. Moreover, to evaluate 
we first need to calculate conditioned on 
fading coefficients vector H and then average the conditional 
BER over H as; 

Pbe{b^^^\b^o\l,d)= f Pbe{bi^^\bi^\l,d,H)f{H)dH, (16) 

JH 


where f{H) is the joint PDE of fading coefficients in H. 
Note that H may be different for each of up-and downlink 
transmissions. We will specify this vector at the end of each 
subsection for each of configurations. 

Different rules can be adopted at the receiver side to decide 
on the transmitted data bits regarding the detected chips. Due 
to the positive nature of incoherent OCDMA systems, the 
presence of interfering users manifests its effect as an increase 
on the photoelectrons count. Therefore, MAI is more likely to 
cause errors on the chips “OEE” rather than on the chips “ON”. 
In order to mitigate this asymmetry in error probabilities 
for chips “ON” and “OEE”, we employ the methodology 
introduced in p5| and p6| . The methodology is based on the 
following decision rule which imposes more restrictions on 
the correct detection of bit “1”: decide that the transmitted bit 
was “1” if all detected chips are “ON” and otherwise recognize 
“0” as the transmitted data bit. In this case, conditional error 
probabilities in ( p^ can be characterized as follows; 


w 

q=l 

W 

ne(0|l,5, H) 

9=1 


(17a) 

(17b) 


in which P^^glg 2 e(l| 0 A: and P^^glg 2 e( 0 |lA, are 

conditional end-to-end chip error rates on the qih transmitted 
chip of the desired user when this chip is “OEE” and “ON”, 
respectively. To obtain approximate end-to-end CERs, we 
make the following assumption: although it is possible to 
detect a chip correctly at the receiver despite of incorrect 
detection in some of the intermediate relays, we neglect these 
lucky events and assume that a chip can be detected correctly 


at the receiver if and only if it is transmitted without any error 
over all intermediate hops. Then, we have; 


AT+l 


= l-n 1, d, 


AT+l 


Piy,2e(0|l,(p,^) = 1-n l-7’ceE(0|l, I, <5, 


i=l 


, (18a) 
, (18b) 


where is the ith hop fading coefficienj^ 

and Pjgl-( 0 | 1 ,are the 
conditional qth chip errors of the ith hop when chips 
“OEE” and “ON” are transmitted, respectively. Therefore, to 
determine ( p^ and ( p^ we begin by calculating CER of each 
intermediate hop. 

1) Conditional CER Analysis for the First Hop: Here we 
calculate the conditional chip error probabilities of the first hop 
using Gaussian approximation, which provides an acceptable 
estimate of the system error rate specially for error probability 
values less than 0.1 | [23| . The first relay’s receiver simply 
integrates over each seconds and compares the result with 
an appropriate threshold to detect the received chip p 6 | . 
In addition to this signal, Poisson distributed dark current 
and background light (with mean photoelectron count rates 
of yd and 75 , respectively) and Gaussian distributed thermal 
noise are added to the received signal count in each 
seconds. Therefore, photo-detected signal after integrate-and- 
dump circuit of the first relay can be modeled as; 


7 = + 


(19) 


in which iq is a vector with W uncorrelated Gaussian dis¬ 
tributed elements each with mean zero and variance of; 


2 _ 2K})TrTc 

“ Re^ 


( 20 ) 


corresponding to the integrated thermal noise p7| . Kij, e, T^, 
and R are Boltzmann’s constant, electron’s charge, the receiver 
equivalent temperature, and load resistance, respectively. On 
the other hand, ^ is an array of Poisson distributed RVs with 
mean ml as; 


Ad) 


(I) 


mi = mi + mi + mi 


Abd) 


( 21 ) 


where based on ( p^ , ml^^^ corresponds to the desired count 
that it is obtained as; 

= ^hybyPcCA^u, (22) 

in which p, h, and / are the photodetector’s quantum effi¬ 
ciency, Planck’s constant, and the optical frequency, respec¬ 
tively. is the aggregated channel loss (due to absorption 
and scattering effects) of the first hop from the desired user to 
the first relay, and u = (1,1,...,!) is a IE-dimensional all-one 
vector. Moreover, mj^^^ that corresponds to the count of the 
undesired interfering users can be expressed as; 


= AfPcPc'^ 

hf 


(I) 


(23) 


^Note that the fading coefficient of the first hop in uplink transmission 
involves more than one coefficient (specifically and as will be 

demonstrated in Section TTT -Al. 







6 


where p is a vector with length W which its q\h element 
is the weighted sum of aq independent log-normal 
RVs, corresponding to the sum of interfering users’ fading 
coefficients, i.e., 


= E ^ 

neAq 


(1)1(1) 


(24) 


in which is the aggregated channel loss of the first hop 
from the nth user to the first relay. Kq specifies the set of aq 
interfering users on the qih chip. Weighted sum in p4l) can be 


{!) 


approximated as a single log-normal RV 
with log-amplitude mean p (/) and variance as |28 


exp ( 24 ^^) 




= !■"( E - 


neAq 


^2 


= 7 I 11 


-i) \ 

(t 

\^Z^neAq -^n,l j 


(25) 


(26) 


where is the variance of log-amplitude factor = 

2ln(/i4*). Furthermore, for = 0 there is not any interfer¬ 
ence on the gth pulsed mark chip and Kq in ( [24| ) specifies 
an empty set and consequently 4^^ = 0. On the other 
hand, + lb)TcU in ( [2T] ) corresponds to the 

dark current and background light counts vector. We can 
summarize the first user’s uplink average received counts of 
the qi\\ chip conditioned on transmitting chips “OFF” and 
“ON” respectively as; 


m 


(uZ,l) 

0,g 


m 


(wZ,l) 

1,9 




Apt 

hf 


+ (id + 'lb)Tc' 


(27a) 

(27b) 


In order to evaluate the first hop’s conditional CER, we 
apply Gaussian approximation which approximates all Poisson 
distributed RVs with Gaussian distributed RVs, where each 
of these Gaussian variables has equal mean and variance. 
The receiver of the first hop selects the threshold based on 
Gaussian approximation as &g = {(TqEi aiEo)/{ai + ao), 

where ao and ai are standard deviations, and Eq and Ei 
are means of the Gaussian distributed decision variable when 
transmitted chips are “OFF” and “ON”, respectively 1^ . 
Note that the first hop’s receiver is synchronized with the 
desired user and hence is imperceptive to the transmitted 
chips of interfering users. Also we assume that the chan¬ 
nel state information (CSI) is available at the receiver, i.e., 
the first hop’s receiver perfectly knows h^ii . Therefore, 
based on Eqs. (p^-([27]) and Gaussian approximation | [23| , 
Eq = (jd + lb)Tc, El = + Eq, (Jq = 

\/Eq + aK and ai = \/Ei + aK . Also ©o simplifies to 
0 , = ^{E,+ai^){E, + ai^) - 
Now we can obtain the first hop gth chip conditional error 


probabilities as; 

Pir')(l|0, 1, a, a, 


= Q 


Qo - m, 


{ul,l) 

0,9 


(uZ,l) 


, (28a) 


pE“')(0|1, l, a, = Pr(J7y> < 0,1?, a, 


= Q 


(ul,l) 

h,q 


Qn 


(ul,l) 

P.9 


, (28b) 


in which and are Gaussian distributed decision 

variables of the first receiver, when the transmitted chips 
are “OEE” and “ON”, respectively. has mean 

and variance has mean 

variance that = EQ-\-Eq^^ and 

^(nu) = El P e[^\ where e[^^ = Therefore, 

MAI imposes similar effect to additive noise. 

2 ) Conditional CER Analysis for the i-th Hop (i=2,...,N-\-l): 
Eor the remaining hops MAI does not affect the chip detection 
process, but underwater turbulence as well as different noises 
do affect the system performance and cause errors on the chip 
detection processes. Similar to ( [27] ) we can express the qth 
chip average received counts of the ith hop conditioned on 
transmitting chips “OEE” and “ON” respectively as; 

^o!q = hd + lb)Tc, (29a) 

+ {id + lb)Tc, (29b) 


where is the ith hop channel loss. Consequently, gth chip 
conditional error probabilities of the ith hop can be determined 
using Gaussian approximation as 


I, a, U^) = (1|0, l, a, ftW) 



(30) 


Note that ^2^ and ( [30] ) can be applied to ( [TS] ) to achieve 
the end-to-end conditional CERs, while insertin g ([T^ in ( [TT] ) 

results into conditional BERs. Moreover, Pbe{bi^^\b^o\ba) 
can be determined by averaging ( [Tt] ) over H as Here, 
for uplink transmission the conditional BER equations in ( [Tt] ) 
can involve maximum of N W 1 fading coefficients as 
H = , I3[^\ I3[P Eor this case 

joint PDE of f{H) in ( [T^ simplifies to f{H) = f{h^ii) x 
Yl^i fiPq^^) ^ independency of 

fading coefficients. In fact the size of H depends on a and 
for each a Eq.® can be evaluated through a J-dimensional 
integral, where J G {N -h 1, + 2,..., N PW P 1}[^ 


^Note that each Ah hop (i = 2, 3,1) introduces a fading coefficient 
of while the first hop presents at least one fading coefficient of h^i . 
Moreover, the first hop introduces a new fading coefficient of for each 
chip mark q = 1,2, as demonstrated in Section III-Al. Depending 

on d, H has a minimum dimension of A/^ + 1 (e.g., when a = (0, 0, 0) 
and all /3g^^s are zero) and maximum dimension of N PW P 1 p.g., when 
d = (1, 2, 3) and each chip experiences a new fading coefficient of 4^^)- 
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B. Downlink BER Analysis 

As it has been shown in | [24| , MAI of the synchronous 
downlink transmission can be eliminated when the number of 
concurrent users satisfy the following condition; 


M < 



(31) 


and transmitting bit bo at the 0 th time slot, and can be 
calculated using Gaussian approximation a^ 





(35) 


In other words, OOC codes would not overlap and their 
internal multiplication, < >, is equal to zero for 

i ^ j. Therefore, qth pulsed mark chip conditional error 
probabilities of the first hop during downlink transmission can 
be determined similar to those of intermediate hops (in the 
previous subsection) as; 




^(1) 


iq,dl). 


;(l)^ 



(32) 


where and are the mean received counts of the ith 
hop when bits “ 0 ” and “ 1 ” are transmitted, and are respectively 
characterized as follows; 

= {id + lb)Tb, (36a) 

+ (7d + 7b)Tb, (36b) 

in which and are respectively the transmitted power 

per bit and channel loss of the ith hop. It can be shown that due 

to the independency of fading coefficients in H, ( [34| ) simplifies 
to; 


where and are the gth chip average received 

counts of the first hop conditioned on transmitting chips “OFF” 
and “ON”, and are respectively defined as; 

= {ld + lb)Tc, (33a) 

+ (id + lb)Tc. (33b) 

Conditional CER for each of the other hops is exactly 
that of uplink transmission in ( [30| ). Eventually, relying on 
( [30| ) and ( [3^ the end-to-end conditional CER and BER 
of downlink transmission can be evaluate d si milar to 
and ([nl., respectively. Eurthermore, P^g ( 6 o^^| 6 q^^) can be 
determined by averaging <[ni> over H as P^e(^o^^ |^o^^) = 
Jfj Pi)e{b^o^\h^Q\H)f{H)dH. It is obvious that for downlink 
transmission H is defined as a {N + 1)-dimensional vector 
as P = and joint PDE of f{H) 

simplifies to f{H) = f{h^ii) x 
independency of fading coefficients. 


C. Single-User Multi-Hop Point-to-Point UWOC Link 

As a special case, the BER performance of a single-user 
relay-assisted point-to-point UWOC link with bit-detect-and 
forward strategy can easily be achieved from those obtained 
in Sections III-A and III-B, by replacing M = 1,P = IU = 1, 
and Tg = T 5 . In this case based on (p^-(p^, the end-to-end 
bit error probability can be calculated as; 


Pbe—p2p{^o\^o) — 

/* / ^+1 _ 

/_ 1 - n U - Pbe-i{bo\bo,h^^^) 

\ i=i 


f{H)dH, (34) 


in which bo is the 0th time slot transmitted bit. H and f{H) 
are defined as H = {h^^\h^'^\ and f{H) = 

respectively. Pbe-i{bo\bo,h^"'’) is the ith hop 
bit error probability, conditioned on fading coefficient of 


AT+l 

Pbe-p2p{Mbo) = 1 - n [1 - Pbe-i(b^\bo)] , (37) 

i=l 

where Pbg_^( 6 o|^o) is the ith hop BER conditioned on trans¬ 
mitting bit bo at the 0 th time slot which can be evaluated as 
Pbe-i{Vo\bo) = J~^^oPbe-i{h\bo,h(^^)f{h^^^)dU^b Eq. (|^ 
indicates that the average BER in this case can be obtained 
through the combination of one-dimensional integrals, while it 
was obtained through J-dimensional (or (P + l)-dimensional) 
integral for uplink (or downlink) transmission as we saw in 
Section III-A (or III-B). 

As we demonstrated the conditional chip/bit error prob¬ 
ability of each hop can be considered as a building block 
for evaluating the end-to-end BERs of each uplink, downlink 
or single-user transmissions. This error probability can be 
obtained using either photon-counting method p3| , p9| or 
general classical methods in digital communications, similar 
to Q, 1^ , 1^ . The second approach approximates the 
optical signal as a constant coefficient and all of the noise 
components with an additive Gaussian distributed RV, while 
photon-counting method, which we used in this paper, is a 
more accurate method since it considers the incoming optical 
signal with Poisson distribution. However, it has been shown 
in p9| that both of the above mentioned methods provide 
relatively the same results. On the other hand, different photon¬ 
counting methods exist, e.g., saddle-point and Gaussian ap¬ 
proximations Saddle-point is a sufficiently accurate 

method, but involves some complex nonlinear equations and 
hence its evaluation may be time-consuming. However, in p9| 
the BER of UWOC systems is obtained using both saddle- 
point and Gaussian approximations and it has been concluded 
that Gaussian approximation provides relatively similar results 
to saddle-point approximation. Therefore, in this paper we 
used Gaussian approximation as a relatively accurate and fast 
photon-counting method. Moreover, similar to Q, pO} we can 

^It is worth mentioning that all Eqs. ( [^ , and can also directly 
be derived from by replacing I = 0, d = (0, = 0 and 

eP=o. 
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TABLE I 

Some of the important parameters for the channel simulation 

AND NOISE CHARACTERIZATION. 


TABLE II 

Channel loss and scintillation index for several link ranges 
WITH 6div = 0.02° AND Tc = 10 ns. 


Coefficient 

Symbol 

Value 

Half angle field of view 

FOV 

40° 

Receiver aperture diameter 

Do 

20 cm 

Source wavelength 

X 

532 nm 

Water refractive index 

n 

1.331 

Photon weight threshold at the 

Wth 

1 

o 

1—1 

receiver 



Quantum efficiency 

r] 

0.8 

Equivalent temperature 

Te 

290 K 

Load resistance 

Rl 

100 n 

Dark current 

Tc 

1.226 X 10“° A 

Background mean count rate 

lb 

1.206 X 10^° 1/s 

Rate of dissipation of mean-square 

XT 

10“^ K^/s 

temperature 



Rate of dissipation of turbulent 

e 

5 X 10-^ 


kinetic energy per unit mass of fluid 


Relative strength of temperature and w 

salinity fluctuations 

-3.5 

Normalized parameters of the (A, ©) 

Gaussian beam 

(0,1) 


use Gauss-Hermite quadrature formula (D to approximate J- 
dimensional integrals with J-dimensional finite series. 

IV. Numerical Results 

In this section, we investigate BER performance of the 
proposed system for different scenarios. To do so, we consider 
an underwater network established in clear ocean. For this 
water type parameters a, b, and c are 0.114 m“^, 0.037 
m“^, and 0.151 m“^, respectively. We assume that all users 
are equidistant to the destination (or to the first relay in the 
relay-assisted configuration). Further, we assume that when 
N intermediate relays are employed, all of the hops have 
the same link length of do/{N 1) m and therefore the 
same channel loss and fading log-amplitude variance. To 
evaluate the channel loss of each hop, we have simulated the 
channel impulse response based on MC approach similar to 
Q, |T^ . Moreover, to obtain fading coefficient of each hop 
and dependency of these coefficients to the link range we 
consider specific values for the water turbulence parameters 
and numerically calculate the scintillation index for various 
link ranges similar to GD- Table I displays some important 
parameters of the system for the channel simulation as well 
as noise characterization. 

As we mentioned in Section II-B, from the practical aspect 
of view it is more preferable to use diffusive links for the first 
and/or last hop(s) to reduce the pointing and tracking diffi¬ 
culties in these hops. In order to include both collimated and 
diffusive links, we first deal with collimated links developed 
by narrow laser beams. Then, at the remaining of this section, 
we investigate the effect of diffusive beams on the system BER 
and also on the validity of our assumption in Eq. GD- Table 
II summarizes the simulation results for the channel loss and 
scintillation index in several link ranges with the source full 
beam divergence angle of Odiv = 0.02^ and the chip duration 
time of Tc = 10 ns. 

In our numerical results we use OOC signature with the 
following parameters, namely (T, IE, Aa,Ac) = (50,3,1,1). 


Link 

range 

Channel loss 
coefficient, L 

Scintillation 
index, aj 

Log-amplitude 

variance, 

90 m 

3.99 X 10“^ 

0.9738 

0.17 

70 m 

7.812 X 10-° 

0.616 

0.12 

45 m 

3.135 X 10“^ 

0.271 

0.06 

30 m 

3.1 X 10“^ 

0.1248 

0.029 

22.5 m 

9.4 X 10-=^ 

0.071 

0.017 

18 m 

18.2 X 10-^ 

0.0452 

0.011 


Table III summarizes the results of the interference char¬ 
acterization for M = 5 simultaneous users, each with the 
aforementioned OOC parameters. First column presents the 
possible number of interferences on the desired user while the 
second column indicates probability of having I interferers, 

defined as Pi{l) = (^1 - . As it 

can be observed, larger number of interferers occur with less 
probability. Third column is the total number of possible 
interference patterns a, when I interference occurs on the 
W chips of the desired user. In other words, this column 
is the number of as satisfying condition aq = /, 

which can be obtained as Fourth column illustrates 

possible interference patterns for each /, while probability of 
each pattern conditioned on having I interferers, P^(a|/) = 
TT./ -i-rw / —G. is determined in the fifth column. Note that 

Uq = l{aqV. 

different users’ interference is more likely to occur on the 
different chips instead of occurring on a certain chip, e.g., 
interference pattern of a = (0,0,3) is less probable than 
a = (1,1,1). Sixth column represents number of similar 
as for each row of the fourth column, i.e., number of as 
that have the same components as the interference pattern 
of the fourth column. For example, we have 5 interference 
patterns similar to (0,1,3), namely (0,3,1), (1,0,3), (3,0,1), 
(1,3,0) and (3,1,0). Each of these similar patterns has the 
same conditional probability of occurrence as the fifth column. 
Moreover, dimension of H in ( p^ for uplink transmission 
is summarized in the seventh column. We should emphasize 
that crosscorrelation constraint of Ac = 1 imposes that each 
interfering user can make interference maximum on a one chip 
of the desired user |[^. For example a = (1,1,1) indicates 
three different interferers, each interfering on a certain chip of 
the desired user; resulting into three new fading coefficients 
l3i\ and 

Fig. 3 illustrates the fading-free impulse response of the 
channel for various link ranges, simulated based on MC 
method. As it can be seen, even for link range of do = 90 m 
the clear ocean channel impulse response can approximately 
be modeled by an ideal delta function. This behaviour strongly 
guarantees GT), particularly for Tc = 10 ns. Therefore, 
we consider the chip duration time as Tc = 10 ns, and 
consequently each user can communicate with data rate of 
Rh = 2 Mbps. 

Fig. 4 indicates the scintillation index cr| and the corre¬ 
sponding log-amplitude variance (7\ of the channel in terms 
of link range, with respect to the turbulence parameters in 
Table I. Furthermore, effects of increase on these parameters 
is investigated in this figure. As it can be seen, increases on 
































9 


TABLE III 

Interference characterization for M = 5 concurrent users 
WITH OOC PARAMETERS OF (F, W, Xa, Ac) = (50, 3, 1, 1). 


1 

m) 

Total # 
of as, 

('t^) 

d 

Ps{S\l) 

# of 
similar 
ds 

Dimension 
of H in 

06 ) 

0 

0.6857 

1 

(0,0,0) 

1 

0 

N + 1 

1 

0.2713 

3 

(0,0,1) 

1/3 

2 

77 + 2 

2 

0.0402 

6 

(0,0,2) 

1/9 

2 

A^ + 2 

(0,1,1) 

2/9 

2 

N + 3 

3 

0.0027 

10 

(0,0,3) 

1/27 

2 

77 + 2 

(0,1,2) 

1/9 

5 

77 + 3 

(1.1,1) 

2/9 

0 

Ar + 4 

4 

6.561 

xlO-5 

15 

(0,0,4) 

1/81 

2 

77 + 2 

(0,1,3) 

4/81 

5 

77 + 3 

(0,2,2) 

2/27 

2 

A^ + 3 

(1,1,2) 

4/27 

2 

A^ + 4 



- 90 m 



4 4.1 4.2 4.3 

time, t [s] X 10"" 


O'— 
1.95 


0.01 

0.005 

0 


- 45 m 


2 2.05 

time, t [s] 



- 22.5 m 



0 ^ 

2.8 


3 3.2 3.4 3.6 

time, t [s] X iQ-^ 


1 1.05 1.1 1.15 

time, t [s] ^ 


^ 0.02 
1 0.01 
0 


- 30 m 


1.3 


1.35 1.4 1.45 1.5 

time, t [s] ^ iQ-^ 



1 18 m 1 



8.1 8.2 8.3 

time, t [s] ^ iQ-s 


Fig. 3. Fading-free impulse response of the clear ocean channel for various 
link ranges, do = 90 m, 70 m, 45 m, 30 m, 22.5 m, and 18 m. 


XT and w results into stronger fading while smaller value for 
the scintillation index is expected for larger 5 . 

In order to perform the numerical simulations we generate 
10 ^ Bernoulli RVs with PDF of Pr(6^"5 = + 

— 1) for each user, where S(.) is Dirac delta function. 
Then, we transmit 10^ Poisson distributed RVs with mean 



Fig. 4. Effect of increasing turbulence parameters e, xt and w on the 
scintillation index af. 



-20 -10 0 10 20 

Average transmitted power per bit [dBm] 


Fig. 5. Effect of increasing the number of simultaneous users M and the 
thermal noise variance on the uplink performance of a do = 70 m clear 
ocean link without any intermediate relay node. 

^PcTcb^j^^ for each chip and multiply these RVs with the 
link loss and with 10^ log-normal RVs corresponding to 
the link fading coefficient. At the receiver side we add the 
resulted sequence with 10^ Gaussian distributed RVs that 
correspond to the receiver thermal noise and with 10^ Poisson 
distributed RVs corresponding to the receiver dark current 
and background noise. Afterward, we compare each of the 
10^ resulted decision variables with an appropriate threshold 
to detect the received chips. The detected chip of each node 
would be transmitted via a Poisson distributed RV with mean 
^PcTcCq \ where Cq^ is the qth detected chip of the ith 
node. Pursuing similar procedures, we can detect the received 
chips at the destination. Then, if all the chips of a time slot 
are “ON” we recognize bit “1” as the transmitted data bit and 
vice versa. Finally, the BER can be evaluated comparing the 
transmitted and the decoded data sequences. 

To see effects of increasing the number of simultaneous 
users on the uplink performance, in Fig. 5 the uplink BER 
of a do = 70 m clear ocean link without any intermediate 
relay is depicted for M = 1, 3, 4 and 5, = 2 Mbps, and 

two different values of the thermal noise variance, namely 
= 10^ and 10^. As it can be seen, when the number 
of simultaneous users is smaller or equal to the code weight 
(M < W) MAI does not affect the system performance. 
However, as the number of users violates the code weight, the 
uplink BER saturates to a constant value. For example, when 
M = 4 the main contribution of BER is due to the interference 
pattern of a = (1,1,1). In fact for this interference pattern 
7^6e(l|0) ^ 0 ^^d Pe ^ = 3, (T = 

(1,1,1)) = i X 1.62 X 10-^ = 8.1 X 10-^ Note that as the 
number of users increases the saturation value increases too. 
Also larger BERs are expected for larger values of the thermal 
noise variance. 

However, the decision rule of Eq. appropriately 

mitigates the bit error probabilities asymmetry for weak 
to moderate SNR regimes, at high SNR regimes the se¬ 
vere increase on the photoelectrons count of the desired 
user’s chips, caused by high interference power of unde¬ 
sired users, considerably increases the difference between 
error probabilities for chips “ON” and “OFF”. In other 
words, at high SNR regimes where the system is MAI- 
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Fig. 6. Uplink BER of an OCDMA-based underwater network without any 
intermediate node. M = 5, = 3.12 x 10^, Ri) = 2 Mbps, and various 

communication ranges do = 90 m, 70 m, 45 m, 30 m, 22.5 m and 18 m. 

limited > ^ie-e 2 e( 0 |l’ 

hence Pi)e{l\0,l,d, H) > P^e(0|l: The aforemen¬ 

tioned asymmetry is mainly due to the simple threshold value 
of the receiver that only takes into account the transmitted 
power of the desired user and its channel condition. In general, 
maximum likelihood (ML) hypothesis should be considered 
to derive the optimum threshold value and hence the op¬ 
timum receiver structure. Undoubtedly, the exact BER of 
the OCDMA system with optimum threshold value (which 
relates to the number of interfering users and their transmitted 
power and channel condition) would be less than the above 
mentioned simple receiver’s BER. Therefore, the receiver 
structure adopted in this paper provides a simple decision rule 
and hence an upper bound on the BER 

Eig. 6 illustrates uplink BER of the underwater network 
for various link ranges and M = 5 simultaneous users that 
transmit OCDMA-based data to a destination directly without 
any intermediate relay. Also based on Eq. ( |2Q| ) and values 
of the parameters in Table I, the thermal noise variance is 
= 3.12 X 10^. As it can be observed, increasing the link 
range severely degrades the performance. This is reasonable 
since the degrading effects of absorption, scattering and fading 
increase rapidly with distance. Hence, the idea behind multi¬ 
hop transmission becomes apparent; divide the long commu¬ 
nication distance to shorter ones (by means of intermediate re¬ 
lays) each with much reduced degrading effects of absorption, 
scattering and fading. Moreover, the well matches between the 
analytical and numerical simulation results show the accuracy 
of our derived expressions. However, a few discrepancy exists 
at low error rates, mainly due to numerical calculation of 
multi-dimensional integrals of analytical expressions. We also 
calculated the uplink BER without approximation to the sum 
of log-normal RVs, using M-dimensional integrals. As it 
can be seen, the well match between the results confirms 
the accuracy of the mentioned approximation in Eqs. 

( [2^ . In general, this approximation is more beneficial when 
M > IU +1, since it converts (M ^A/")-dimensional integrals 
to (lU -f 1 -b A^)-dimensional counterparts. 

Eig. 7 depicts uplink BER of the relay-assisted UWOC 
network where M = 5 OCDMA-based simultaneous users 
are communicating with a destination located 90 m away 



Fig. 7. Uplink BER of the relay-assisted OCDMA-based underwater network 
with M = 5 users and end-to-end distance of do = 90 m. Rj, =2 Mbps, 
0-2, = 3.12 X 10^ and AT = 0, 1, 2, 3 and 4. 



Fig. 8. Downlink BER of the relay-assisted OCDMA-based underwater 
network with M = 5 simultaneous users, do = 90 m, = 2 Mbps, 
= 3.12 X 10'^, and AT = 0, 1, 2, 3 and 4. 

in clear ocean. In this scenario, we assume that all of the 
transmitters transmit their corresponding chips with the same 
multi-hop chip power of where Pb^avg 

is the average transmitted power per bit of each user. As 
it is shown in Eig. 7, significant performance improvement 
can be achieved by utilizing intermediate relays, e.g., 30 dB 
improvement in the BER of 10“^ by employing dual-hop 
transmission. Comparing Eigs. 6 and 7 shows that in multi¬ 
hop transmission impairment of each hop is a bottleneck on 
the overall performance. In other words, a multi-hop link 
with N intermediate relays and do m end-to-end distance 
can approximately achieve the performance of a single-hop 
link with do/(A^ + 1) m communication distance and by 
approximately 10 logio(^ + 1) larger average transmitted 
power per bit (here do = 90 m and A^ = 0,1,..., 4). 

Eig. 8 indicates downlink BER of the relay-assisted 
OCDMA-based UWOC network. Eor this case OBTS trans¬ 
mits synchronous data to a collection of M = 5 users (located 
90 m away in clear ocean), which satisfies the inequality 
in (ju). Therefore, downlink transmission does not suffer 
from MAI. Consequently, downlink BER can surpass that of 
uplink (when M > lU + 1), as it is apparent by comparing 
Eigs. 7 and 8. Moreover, similar to uplink transmission multi¬ 
hop downlink transmission can yield significant performance 
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Fig. 9. BER performance of a single-user relay-assisted point-to-point 
UWOC system, do = 90 m, Rb = 100 Mbps, = 3.12 x lO'^, N = 0, 
1, 2, 3 and 4. 


X 10"‘' 



time, t [s] 



time, t [s] 


xl0“® 





Fig. 10. Beneficial application of the multi-hop transmission (serial relaying) 
in expanding the viable communication range; BER of a single-user relay- 
assisted point-to-point UWOC system with N = 3 intermediate relays, in 
several link ranges, i.e., do = 280 m, 180 m, 120 m, 90 m and 72 m. 

improvement, e.g., 32.5 dB in the BER of 10“^ by means of 
only one intermediate relay. 

Fig. 9 shows BER performance of a single-user relay- 
assisted point-to point UWOC system with data rate of 100 
Mbps in a 90 m clear ocean link. As it can be seen multi-hop 
transmission significantly improves the system performance, 
e.g., 32 dB improvement in the BER of 10“^ is achieved 
using a dual-hop transmission. We should emphasize that the 
negligible discrepancy between the simulation and analytical 
results of relay-assisted topology at high BER regimes is 
caused by the assumption in ( p^ , which provides a relatively 
tight upper bound on the system BER. 

In order to more clearly illustrate the beneficial application 
of multi-hop transmission in extending the viable communi¬ 
cation distance, we consider a single-user relay-assisted point- 
to-point UWOC system with N = 3 intermediate relays. 
Fig. 10 depicts the BER of such a system in several link 
ranges. As it can bee seen, by means of multi-hop transmission 
we can attain comparatively large communication distances 
with realistic transmit powers and acceptable BERs. For 
instance, we can communicate through a 280 m clear ocean 
link with an appropriate performance and a realistic average 
transmitted power per bit. These achievable communication 
distances substantially differ from the reported distances in the 


Eig. 11. Fading-free impulse response of the clear ocean link with transmitter 
beam divergence angle of Odiv = 5° and various link lengths, namely do = 
90 m, 45 m, 30 m and 20 m. 

existing literature (i.e., ‘'typically shorter than 100 m”). Our 
numerical results confirm that employment of intermediate 
relay nodes is of utmost importance for underwater optical 
communication. Appropriate number of relays can be obtained 
using the network designing policies. For example if the 
system is power constrained, we can first obtain the BER of the 
system under a set of conditions versus the transmitted power 
and then we can choose the number of relays such that an 
acceptable performance achieves using a predetermined total 
transmitted power. Moreover, the largest link span depends 
on the system parameters and for a given set of conditions 
and system parameters it can easily be determined using our 
derived expressions for the system BER. 

However, our previous numerical results dealt with colli¬ 
mated laser beams for all hops, at the remaining of this section 
we investigate the effect of diffusive links. Fig. 11 illustrates 
the MC simulation results for the fading-free impulse response 
of UWOC channel with different link ranges, namely do = 90 
m, 45 m, 30 m and 20 m. Also the transmitter full beam 
divergence angle is assumed to be Odiv = 5^. As it can be seen, 
even for 90 m clear ocean link the channel spreading time is 
negligible with respect to the chip duration time. Based on 
( p^ , the channel loss coefficients for to the above mentioned 
link lengths are 5.9224 x 10“^ 1.8575 x 10-^ 2.5194 x 10“^ 
and 1.6167 X 10 respectively. 

Fig. 12 shows the uplink BER of a relay-assisted OCDMA- 
UWOC network with M = 3 users, N = 2 intermediate 
relays and do = 90 m end-to-end communication distance. In 
this figure we consider diffusive transmitters with Odiv = 5^ 
for the first and last hops. The transmitter of the second hop 
(between two relays) is yet assumed to be collimated with 
Odiv = 0.02^. We also investigate the effect of changing the 
position of intermediate relays on the system BER. for this 
purpose we define do^j as the length of the hop between 
the nodes ith and jth. Here, do,i 2 + do ,23 + do ,34 = 90 
m. We define different modes of operation for the assumed 
network. Mode I specifies an OCDMA-UWOC network with 
equal hop lengths as do, 12 = do,23 = do,34 = 30 m, same 
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Fig. 12. Uplink BER of a relay-assisted OCDMA-UWOC network with 
M = 3, A/' = 2, do = 90 m and = 3.12 x 10^. The effects of diffusive 
links for the first and last hops as well as changing the position of intermediate 
relays are investigated. 

TABLE IV 

Channel loss and scintillation index lor several link ranges 
CORRESPONDING TO ElG. 12. 


Link 

range 

Full beam 
divergence 
angle, 

^div 

Channel loss 
coefficient, L 

Scintillation 
index, aj 

Log- 

amplitude 

variance. 

30 m 


2.5194x10“^ 

0.1248 

0.0294 

35 m 

0.02*' 

1.4722x10-^ 

0.1681 

0.03884 

27.5 m 

5^ 

3.8713x10-^ 

0.1054 

0.02505 

40 m 

0 . 02 'j 

6.9336 xl0-"‘ 

0.2169 

0.04907 

25 m 

5"^ 

6.1802x10-^ 

0.0873 

0.02093 

22.5 m 

5^ 

1.0015x10-^ 

0.0712 

0.0172 

50 m 

0 . 02 'j 

1.5282x10-"' 

0.3303 

0.07136 

20 m 

5"^ 

1.6167x10-'* 

0.0559 

0.01361 


FOV = 40^ for all of the receivers and equal divergence angle 
of Odiv = 0.02^ for all of the transmitters. Mode II relates to 
the same parameters as Mode I and only the transmitters of 
the first and last hops are assumed diffusive with Odiv = 5^. 
For Modes III, IV, V and VI we assume exactly the same 
parameters as Mode II, except that the effect of changing 
the position of intermediate nodes is investigated. Specifically, 
6 ^ 0,12 = <^0,34 = 27.5 m and do, 23 = 35 m for Mode III, 
<^0,12 = <^0,34 = 25 m and do, 23 = 40 m for Mode IV, 
do,12 = <^0,34 = 22.5 m and do,23 = 45 m for Mode V, 
and do,12 = do ,34 = 20 m and do,23 = 50 m for Mode VI. 
The channel loss coefficient as well as fading statistics for 
different link rages and various transmitter beam divergence 
angles are calculated and summarized in Table IV. As it can be 
seen, increasing the beam divergence angle of the first and last 
hops’ transmitters degrades the system performance; however 
it is important from the practical point of view. Moreover, 
our numerical results show that we can achieve a better 
performance by changing the position of intermediate relays, 
i.e., there is an optimal place for intermediate relays from the 
BER viewpoint. For example. Mode IV has approximately 
7.5 dB better performance than Mode VI at the BER of 
10“^. It is worth noting that as the total link range and the 
beam divergence angles of diffusive transmitters increase, the 
differences between performances of various modes increase 
too and hence the optimal relay placement achieves more 
importance. 


V. Conclusions 

In this paper, we studied the BER performance of optical 
code division multiple access (OCDMA)-based underwater 
wireless networks for both up-and downlink transmissions. 
We considered all the disturbing effects of the medium to 
better model the underwater channel. To benefit from the 
incremental dependency of absorption, scattering and fading to 
the underwater link range, we considered relay-assisted trans¬ 
mission to divide a relatively long communication distance to 
shorter ones, each with much reduced absorption, scattering 
and fading effects. For the sake of simplicity, we assumed 
chip detect-and-forward strategy at the relay nodes. In order 
to evaluate the chip error probability of each hop and then the 
end-to-end BERs, we applied Gaussian approximation which 
is based on photon-counting methods. We found that when 
the number of simultaneous users is smaller or equal to the 
code weight (M < W), multiple access interference (MAI) 
does not affect the system performance and therefore both up- 
and downlink transmissions have the same end-to-end BERs. 
However, as the number of users violates the code weight 
uplink BER saturates to a constant value; whereas downlink 
BER does not change, since MAI is removed by synchronous 
downlink transmission. Furthermore, BER of a single-user 
relay-assisted point-to-point UWOC system obtained as a 
special case of the proposed relay-assisted OCDMA network. 
Based on our numerical results multi-hop transmission con¬ 
siderably improves the performance, especially for longer link 
ranges, e.g., 32.5 dB performance enhancement obtained via 
only one intermediate relay node for downlink transmission in 
a do = 90 m clear ocean link and at the BER of 10“^. Also we 
numerically showed that there could be an optimal place for 
relay nodes from the BER point of view. Moreover, our sim¬ 
ulation results confirmed that the employment of intermediate 
relay nodes is of utmost importance for underwater wireless 
optical communication, especially for larger distances. 
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